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ABSTRACT 

We report the results of the first sensitive L-hand survey of the intermediate age 
(2.5 - 30 Myr) clusters NGC 2264, NGC 2362 and NGC 1960. We use JHKL colors 
to obtain a census of the circumstellar disk fractions in each cluster. We find disk 
fractions of 52% ± 10%, 12% ± 4% and 3% ± 3% for the three clusters respectively. 
Together with our previously published JHKL investigations of the younger NGC 
2024, Trapezium, and IC 348 clusters, we have completed the first systematic and 
homogenous survey for circumstellar disks in a sample of young clusters that both span 
a significant range in age (0.3 - 30 Myr) and contain statistically significant numbers 
of stars whose masses span nearly the entire stellar mass spectrum. Analysis of the 
combined survey indicates that the cluster disk fraction is initially very high (> 80%) 
and rapidly decreases with increasing cluster age, such that half the stars within the 
clusters lose their disks in ^ 3 Myr. Moreover, these observations yield an overall disk 
lifetime of ~ 6 Myr in the surveyed cluster sample. This is the timescale for essentially 
all the stars in a cluster to lose their disks. This should set a meaningful constraint for 
the planet building timescale in stellar clusters. The implications of these results for 
current theories of planet formation are briefly discussed. 
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1. Introduction 

The fraction of stars presently being formed in the Galaxy that will develop accompanying 
planetary systems, ips, is the product of the fraction of stars initially surrounded by disks, fj, and 
the fraction of original disks that will evolve to form planets, fp, i.e., ips = f « x fp. Both quantities 
(i.e., fj and fp) may be functions of stellar mass, among other astrophysical parameters. The 
initial disk fraction (fj) around newly formed stars can be determined directly from observations. 
However, to determine the probability of planet formation (fp) in such disks requires both empirical 
and theoretical knowledge of the planet formation process. As a first step, empirical knowledge 
of disk evolution can be derived from observations of young stars of varying age and used to 
determine disk lifetimes and thus the amount of time available for planet building around these 
objects. Such information should provide critical tests of planet formation theories and ultimately 
constraints for fp. 

Infrared studies of individual star forming regions such as Ophiuchus, Taurus, Mon R2, 
NGC 1333, Orion, IC 348, Sco-Cen and NGC 2362 have long suggested both that the initial disk 
frequency is relatively high (fj > 50%; Wilking, Lada &: Young 1989 , ^tauffer et al. 1994 , Kenyon 



fc Hartmann 1995| , Lada, Alves fc Lada 1996 , Carpenter et al. 1997 ) and that disk lifetimes are 



relatively short, (3-15 x 10^ yrs.; Strom et al. 198^ , Skrutskie et al. 1990| , ^trom, Edwards, fc 



Skrutskie 1993| , [Lada &: Lada 1995| , [Brandner et al. 2000| , [Alves et al. 2001| ). However, these 



studies suffer from a number of limitations, such as incomplete samples of stars, inhomogenous 
data, uncertain age determinations, etc., which make it difficult to draw firm conclusions about 
either the exact value of the initial disk frequency or the lifetime of the circumstellar disk phase 
of early stellar evolution. In particular, most of these studies are based on statistics derived from 
near-infrared JHK (1.25, 1.65, 2.2 /im) observations which are not long enough in wavelength to 
produce a complete or unambiguous census of disk bearing stars in a population of young stellar 
objects. Thus, existing studies provide tantilizing, but not particularly strong, constraints to guide 
planet formation theories. 

It is clearly desirable to extend and improve upon existing work and place the question of disk 
evolution on a more firm footing. To do so requires a systematic observational campaign which 
both: 1) is conducted at infrared wavelengths which are sufficiently long to identify disk bearing 
stars with the minimum of ambiguity and 2) includes a statistically meaningful sample of stars 
covering a sufficient range of (relatively well determined) ages. 



2. The L-band Cluster Survey 

In order to determine a more robust estimate of disk fractions and lifetimes in star-forming 
environments, we have embarked on a program to produce a homogenous and systematic L-band 
(3.4 /im) imaging survey of a sample of nearby, young, mostly embedded, clusters. Not only are 
such clusters believed to be the birth sites of most galactic field stars ( [Lada 1992| , Carpenter 2000), 
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but they also individually and collectively contain statistically significant numbers (hundreds) of 
stars, providing a meaningful sampling of essentially the entire stellar mass function. Moreover, 
a survey of clusters enables a more meaningful determination of the disk evolution timescale 
than is possible from a survey of a single star-forming region. This is because: 1) a cluster can 
be characterized by the mean age of all its members, which statistically is a relatively reliable 
indicator of age, and 2) a sample of clusters can be observed which span a significantly wider range 
of age than is typical of stars formed in any individual star- forming region. Finally, our previous 
studies have demonstrated that L-band observations are capable of detecting and identifying 
essentially all the disk bearing stars in a population of young stellar objects ( [Haisch, Lada, & 



Lada 2000; hereafter HLLOO, Lada et al. 2000| , Haisch et al. 2001). The fraction of stars with 



circumstellar disks is given directly by the fraction of stars with measured L-band excess emission. 

Our survey sample consists of six nearby and relatively rich young clusters whose published 
mean ages span a range between 0.3 - 30 Myr. In previous contributions, we have presented the 
results of L-band observations of the NGC 2024, Trapezium, and IC 348 clusters (mean ages of ~ 
0.3 - 2 Myr) which clearly demonstrated that the initial disk frequency is quite high (fj> 80%) 
and independent of stellar mass, at least in these very young clusters (HLLOO, [Lada et al. 2000| , 



Haisch, Lada, k. Lada 200l| ; hereafter HLLOl). These results also indicated that as the age of the 
clusters increased, the disk fraction decreased. The rate of this decrease also appeared to be a 
function of stellar mass, where the disks surrounding high mass stars have shorter lifetimes than 
the disks around stars of lower mass (HLLOl). 

In this paper, we present the results of our L-band imaging surveys of the older NGC 2264, 
NGC 2362, and NGC 1960 clusters (mean ages 2.5 - 30 Myr). We use JHKL colors to obtain a 
census of the infrared excess/circumstellar disk fractions in each cluster. In conjunction with our 
previously published JHKL data for the younger NGC 2024, Trapezium, and IC 348 clusters, we 
are able to investigate the evolution of the cluster disk fraction as a function of age. This study 
represents the first homogenous and systematic i-band survey of a sample of young clusters which 
span an adequate age range to enable a straightforward determination of the timescale for disk 
evolution in clusters. 



3. Observations 

We obtained JHK and L-band observations for each cluster in our study. The L-band data 
for each cluster, with the exception of NGC 2362, were taken with STELIRCAM (fTollestrup fc| 
Willner 1998| ) on the Fred Lawrence Whipple Observatory (FLWO) 1.2 meter telescope on Mt. 
Hopkins, Arizona. The details of the STELIRCAM instrument are discussed elsewhere (HLLOO, 
Lada et al. 2000| ). The basic details of our L-band image reduction procedures for NGC 2024, the 
Trapezium, and IC 348 are described elsewhere (HLLOO, Lada et al 2000, HLLOl) and, except as 
noted below, are the same for the NGC 2264 and NGC 1960 clusters. The L-band data for NGC 
2362 were obtained using NSFCAM at the NASA Infrared Telescope Facility (IRTF) on Mauna 
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Kea, Hawaii (the NSFCAM instrument and image reduction procedures are discussed in HLLOO). 

For the NGC 2264 and NGC 1960 clusters, nine fields were observed in a 3x3 square grid 
which covered an area of ~42 arcmin^. Note that our observations of NGC 2264 include only a 
small part of the entire region (that centered on the southernmost cluster) , while we have surveyed 
most of NGC 1960. The total integration time for each field in these clusters is 18 minutes. In 
NGC 2362, twenty-five fields arranged in a 5x5 square covering an area of ~28 arcmin^ were 
observed (covering most of the cluster). The total integration time for each field was ~ 4.5 
minutes. Within the respective survey regions of each cluster, we detected 56, 45 and 93 sources 
at L band to the completeness limits of our surveys {L = 12.0 in NGC 2264 and NGC 1960; L = 
14.0 in NGC 2362). We find a field star contamination within similar areas to be 6, 8 and 18 stars 
in NGC 2264, NGC 1960 and NGC 2362 respectively. 

The JHK observations of the clusters in our sample were obtained at various observatories. 
The NGC 2264 and Trapezium observations were obtained using STELIRCAM on the 1.2 meter 
telescope at FLWO and are discussed in Muench et al. (2001) and Lada et al. (2000). J, H and 
Kg (2.1 fj,m) observations of the NGC 1960 cluster were taken from the 2MASS data archive 
(http://pegasus.astro.umass.edu). JHK observations of the NGC 2362 cluster were obtained at 



the European Southern Observatory (ESO) 2.3 meter telescope and are discussed in Alves et al. 
(2001). The JHK data for the NGC 2024, and IC 348 clusters were obtained with the NOAO 1.3 
meter telescope and are discussed in HLLOO, and Lada & Lada (1995). 



4. Results 

We have used JHKL color-color diagrams to derive the corresponding background corrected 
infrared excess fractions for the NGC 2264, NGC 2362 and NGC 1960 clusters. An extensive 
discussion of the method we use to determine cluster infrared excess fractions from color-color 
diagrams is given in HLLOO. We find i-band excess fractions of 52 it 10%, 12 it 4%, and 3 it 3% 
for the three clusters respectively. These fractions are significantly lower than those previously 



derived for the younger Trapezium and NGC 2024 clusters (HLLOO, Lada et al. 2000) using similar 
techniques. In Table |l|, we list the derived infrared excess fractions for the entire sample of clusters 
which we have surveyed in JHKL bands. Also listed in the table is the estimated stellar mass at 
the respective completeness limits in each cluster as well as the mean age of each cluster taken 
from the literature with the corresponding reference. The ages of the Trapezium, IC 348 and NGC 
2264 were obtained using PMS tracks of Palla & Stabler (2000). The mean age for NGC 2024 
was adopted from Meyer (1996) and was obtained using the PMS tracks of D'Antona & Mazzitelli 
(1994). On the other hand, the ages of NGC 2362 and NGC 1960 were determined from post-main 
sequence isochrone fitting in the HR diagram. 

In Figure |^ we plot the JHKL excess (disk) fraction as a function of mean age for all the 
clusters in our survey except NGC 2362, where we plot the age based on post-main sequence 
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isochrone fitting as discussed earlier. For comparison we also plot excess fractions in Taurus and 
Chamaeleon I (open triangles), derived from similar JHKL observations in the literature (i.e., 
Kenyon fc Hartmann 1995| , [Kenyon fc Gomez 2001| ). The ages for Taurus and Chamaeleon I 



were obtained from Palla & Stahler (2000). NGC 1960 is not included in the figure since our 
observations of this cluster only extend to > 1 Mq stars, whereas in the other clusters we are 
complete to < 1.0 M©. 

The dot-dashed line in Figure 1 represents a least squares fit to the data obtained in our 
L-band survey (filled triangles). Vertical error bars represent the statistical ^/N errors in our 
derived excess/disk fractions. Horizontal error bars show representative errors of our adopted 
ages. The error bars for the ages of the Trapezium, Taurus, IC 348, Chamaeleon I and NGC 2264 
represent the error in the mean of the individual source ages derived from a single set of PMS 
tracks. In order to estimate the overall systematic uncertainty introduced in using different PMS 
tracks, we calculated the mean age and the standard deviation of the mean age for NGC 2264 (2.6 
± 1.2 Myr) from five different PMS models (|Park et al. 2000| , Palla fc Stahler 200q ). This latter 



quantity illustrates the likely systematic uncertainty introduced by the overall uncertainties in the 
PMS models. This is plotted in Figure 1. For stars with M < 1 Mq and ages < 5 Myr, the overall 
uncertainty in the ages for all regions is likely within about 1 - 1.2 Myr. The plotted error for 
NGC 2024 reflects this uncertainty. The age error for NGC 2362 was adopted from the literature 
QBalona fc Laney 1996| ). 



5. Discussion 

We have completed the first sensitive L-band survey of a sample of young clusters that span a 
sufficent range in age (0.3 - 30 Myr) to enable a meaningful determination of the timescale for disk 
evolution within them. Clusters appear to be characterized by a very high initial disk frequency 
(> 80%) which then sharply decreases with cluster age. Half the disks in a cluster population are 
lost in only about 3 Myr, and the timescale for essentially all the stars to lose their disks appears 
to be about 6 Myrs. 

The precise value of this latter timescale to some extent depends on the derived parameters 
for the NGC 2362 cluster. Our quoted timescale of 6 Myrs could be somewhat of a lower limit 
for two reasons. First, it is possible that a slightly higher disk fraction for NGC 2362 could be 
obtained with deeper L-band observations which better sample the cluster population below 1 Mq. 
Our earlier observations of IC 348 and the Trapezium cluster show that the disk lifetime appears 
to be a function of stellar mass (HLLOl), with higher mass stars losing their disks faster than 
lower mass stars. However, we note that much deeper Ji^i^ observations ( Alves et al. 200l| ) which 



sample the cluster membership down to the hydrogen burning limit yield a JHK disk fraction of 
essentially 0%, giving us confidence in the very low disk fraction derived from our present L-band 
observations. Second, the age of NGC 2362 is dependent on the turnoff age assigned to only one 
star, the O star, r CMa. This star is a multiple system and its luminosity assignment on the HR 
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diagram somewhat uncertain (van Leeuwen & van Genderen 1997). Correction for multiplicity 



would lead to a slightly older age. However, the quoted 1 Myr error in its age likely reflects the 



magnitude of this uncertainty (Balona & Laney 1996). On the other hand, if, for example, the 



errors were twice as large as quoted, the cluster could have an age between 3 Myr and 7 Myr. The 
corresponding age and the overall disk lifetime derived from a least squares fit to the data would 
be between 4 Myr and 8 Myr. Even if the timescale for all disks to be lost was as large as 8 Myr, 
our survey data would still require that half the stars lose their disks on a timescale < 4 Myr. 
Finally, an even older age for NGC 2362 would likely indicate that the decrease in disk fraction 
with time does not follow a single linear fit, that is, after a rapid decline during which most stars 
lost their disks, the disk fraction in clusters would decrease more slowly, with a small number of 
stars (~ 10%) retaining their disks for times comparable to the cluster age. On the other hand, 
we estimate the dynamical age of the S310 HII region which surrounds, and is excited by, r CMa 
to be ~ thii I ^exp ~ 5 X 10^ yrs for rj^// = 50 pc and Vg^-p = 10 km s~^ (e.g. Lada &: Reid 1978 , 



Jonas, Baart, fc Nicolson 1998| ). This is consistent with the turn-off age of the cluster derived 



from the HR diagram, and supports our estimate of ~ 6 Myr for the overall disk lifetime. 

We point out that our L-band observations directly measure the excess caused by the presence 
of small (micron-sized), hot (~ 900 K) dust grains in the inner regions of the circumstellar disk and 
these observations are sensitive to very small amounts (~ lO'^'' gm) of dust. We expect that if there 
is gas in the disk, turbulent motions will always keep significant amounts of small dust particles 
mixed with the gas ( Ruden 1999| ), thus dust should remain a good tracer of gas in the disks as 



they evolve to form planets. Indeed, recent observations of H2 in older debris disks appear to 
confirm this assertion ( [Thi et al. 2001| ). Consequently, stars which did not show infrared excesses 
are likely to be significantly devoid of gas as well as dust. There is also evidence that the presence 
of dust in the inner disk regions is linked with the presence of dust in the outer disk regions (i.e., r 
> 1 AU) where most planet formation is likely to occur. Earlier, HLLOl noted a strong correlation 
between the presence of an L-band excess and millimeter-wave continuum emission (which traces 
the outer disk regions) for stars in the Taurus clouds. Thus the disk evolution timescales derived 
from our L-band survey likely characterize the bulk of the disk material around a young star and, 
consequently, should place important constraints on the timescale allowed for building gas giant 
planets around such stars. 

For example, according to the accepted model of giant planet formation, the core accretion 
model, such planets are built in the following sequence: a high density (rocky) core with a mass 
of ~ 10 Earth masses forms, followed by the accumulation of substantial amounts of nebular 
gas. Typical timescales derived from these models range from a few million years to greater than 
10 Myr depending on the mass accretion rate of solid material, the disk surface density of solid 
material and atmospheric dust abundance (Bodenheimer, Hubickyj, & Lissauer 2000, Lissauei 



2001). A disk lifetime of ~ 6 Myr sets strong constraints on the range of possible values for these 
parameters. The constraints are even more stringent for the highest mass stars, since these stars 
may lose their disks on even shorter timescales (< 2 - 3 Myr; HLLOl). If there is sufficient gas 



-7- 



remaining in the disk after the planet forms, the planet may then migrate inward toward the 
central star ( Lin, Bodenheimer, fc Richardson 1996| ). The discovery of numerous extrasolar planets 



in very close proximity to their central stars (e.g., |Marcy fc Butler 1999| ) would therefore seem 
to require an even more extended disk lifetime. It may be very difficult to both build and then 
migrate planets on such short timescales with standard core accretion models. A competing model 
can produce giant planets on a much shorter (~ 10^ yr) timescale via gravitational instability 
within a relatively massive (Mdisk ~ 0.1 Mq within a radius of 20 AU) protoplanetary disk (Boss 
1998, 2000). The disk lifetimes we measure would be sufficient for giant planets to form via 
gravitational instability. However, whether such models are really viable requires more theoretical 
investigation ([Boss 2000|) . If not, forming giant planets in disks may be more difficult than 
currently envisioned, and the probabilty (fp) of producing giant planets in typical circumstellar 
disks may be relatively small. 

However, we cannot firmly rule out the possibility that the lifetime for the outer disk regions, 
which contain the bulk of the mass available to form planets, may be different from that for the 
inner disk. If the outer disk lifetime is longer than that which we measure for the inner disk, 
models which predict a longer giant planet formation timescale could be accomodated ( p?hi et aL| 



2001 , Lissauer 2001 ). Longer wavelength millimeter and Space Infrared Telescope Facility (SIRTF) 
observations of statistically significant samples of YSOs in a number of young clusters would be 
very beneficial to providing further insight into this question and improving the assessment of the 
overall likelihood of planet formation in circumstellar disks. 
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Fig. 1. — The JHKL excess/disk fraction as a function of mean cluster age. Vertical error bars 
represent the statistical \fN errors in our derived excess/disk fractions. Horizontal error bars show 
representative errors of our adopted ages. See text for discussion. 



Table 1. Background Corrected JHKL Cluster Infrared Excess Fractions 



Cluster JHKL Excess Fraction^ Mass (M©)^ Age (Myr) & Reference^ 



NGC 2024 


112/131 (85% ± 8%) 


0.13 


Trapezium 


312/391 (80% ± 5%) 


0.16 


IC 348 


59/91 (65% ± 8%) 


0.19 


NGC 2264 


26/50 (52% ± 10%) 


0.85 


NGC 2362 


9/75 (12% ± 4%) 


1.00 


NGC 1960 


1/37 (3% ± 3%)4 


1.30 



0.3; iMeyer 1996 



1.5; Palla & Stabler 2000 



2.3; Palla & Stabler 2000 



3.2; Palla & Stabler 2000 



5.0; Balona &: Laney 1996| 



30.0; Barkbatova et al. 1985 



Ks filter was used for tbe observations of NGC 1960. 
^Mass of a star at the completeness limit for eacb cluster. 
^Reference from wbicb cluster age was taken. 

^Tbe source wbicb sbowed an L-band excess is an outflow source superimposed 
on the cluster field and likely not associated with tbe cluster itself (|Magnier et al. 
19961) . 
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